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the behavior of matter at supernuclear densities, e.g., hyperonic degrees 
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itself are tested against the body of observed cooling data. 
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1 Introduction 



A forefront area of research, both experimental and theoretical, concerns the 
exploration of the properties of matter under extreme conditions of tempera- 
ture and/or density and the determination of the equation of state (pressure 
versus density) associated with it. Its knowledge is of key importance for our 
understanding of the physics of the early universe, its evolution to the present 
day, compact stars, various astrophysical phenomena, and laboratory physics 
(for an overview, see, for example, []]]). On the earth, relativistic heavy-ion 
colliders provide the only tool by means of which such matter can be created 
and its properties studied. On the other hand, however, it is well known that 
nature has created a large number of massive stellar objects, i.e., white dwarfs 
and neutron stars, which contain matter in one of the densest forms found in 
the universe. Neutron stars are associated with two classes of astrophysical 
objects - pulsars and compact X-ray sources. Matter in their cores possess 
densities ranging from a few times the density of normal nuclear matter to 
about an order of magnitude higher, depending on star mass. To the present 
day, about 600 pulsars are known, and the discovery rate of new ones is rather 
high. This is accompanied by an impressive growth rate of the body of observed 
pulsar data, like pulsar temperatures determined by the X-ray observatories 
Einstein, EXOSAT, and ROSAT @-@|. 

In this paper, we shall apply a broad collection of modern, field-theoretical 
equations of state (EOS) to the study of the cooling behavior of both neutron 
stars and their strange counterparts - strange matter stars - which should 
exist if 3-flavor strange quark matter is more stable than confined hadronic 
matter. This collection of EOSs was derived under numerous model assump- 
tions about the behavior of superdense stellar matter. To mention several are: 
the many-body technique used to determine the equation of state; the model 
for the nucleon-nucleon interaction; description of electrically charge neutral 
neutron star matter in terms of either only neutrons, neutrons and protons in 
generalized chemical equilibrium {(3 equilibrium) with electrons and muons, 
or nucleons, hyperons and more massive baryon states in f3 equilibrium with 
leptons; behavior of the hyperon coupling strengths in matter, inclusion of 
meson (ir, K) condensation; treatment of the transition of confined hadronic 
matter into quark matter; and assumptions about the true ground state of 
strongly interacting matter (i.e., absolute stability of strange quark matter 
relative to baryon matter). 

The paper is organized as follows. In section |2] we introduce the set of equations 
that govern the cooling behavior of massive stars. The collection of equations 
of state for neutron stars is discussed in section |3|. The physics of strange 
stars and their associated EOS is explained in section |j. The phenomenon 
of superfluidity and the various neutrino emission processes are outlined in 
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sections |] and Our results and conclusions are presented in sections [F] and 
respectively. We summarize the paper in section ^. 



2 Structure equations 

The structure of compact, massive stars is determined by Einstein's field equa- 
tions, 

= 8ttT^ , (1) 

where G^ u denotes the Einstein tensor and T^ v is the energy-momentum ten- 
sor. For spherically symmetric, non-rotating stellar configurations the line el- 
ement is given by (Schwarzschild metric) 

ds 2 = _ ^ dt 2 + e 2A dr 2 + r 2 dQ 2 + r 2 ^2 q ^2 ? (2) 



where <f>(r) and A(r) denote the radially dependent metric functions. 

Equations (|l|) and (|2]) combined with the law of energy-momentum conserva- 
tion, T^ u = 0, lead to 



*L = I ( 3 ) 

dN 4vrr%e A ' 1 ; 

dM p _ A 



(4) 
(5) 



dN n 

d(f) _G (§ r 3 P + M) e 



dN Anc 2 r 4 n 
dP 



which describe the structure of static stars that are in general-relativistic hy- 
drostatic equilibrium (Tolman-Oppenheimer-Volkoff equations, cf. |5|,§]). The 
metric function A is determined by 



2GM )~ l . (7) 



rc 2 



The baryon number N is introduced as an integration variable, since this is 
the only conserved quantity in relativistic field theory. M(N) is the total mass 
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within a sphere that contains N baryons. Its surface area is equal to 47rr 2 (iV). 
The quantity P{N) denotes pressure as a function of baryon number, p(P(N)) 
is the total energy density, and n(P(N)) the baryon number density. Finally, 
we mention that the functions p{P) and n(P) are related to each other by the 
equation of state (EOS). 

The corresponding boundary conditions are given by 



r (iV = 0) = 0, (8) 
M(iV = 0) = 0, (9) 

2GM(N = N m y 
r(N = N m )c 2 

p(P(N = N m )) = p m . (11) 



«" = Ar m ) = ll» 1"=^ . do) 



The quantity N m is defined by the relation 

M(N = N m ) = M m , (12) 



where M m denotes a given star mass. 

The photosphere of the star is treated separately from the structure equations 
(@) _ @ f° r t wo reasons: Firstly, the equation of state of the photosphere de- 
pends much more sensitively on temperature than the core's equation of state, 
which can be treated at zero temperature; secondly, because of a scaling be- 
havior (see, e.g., J7]|§) the structure equations of the photosphere need to be 
solved only once, for a properly chosen star mass. The interface between crust 
and photosphere is located at the density p m . We shall use the value given by 
Van Riper [g], p m = 10 10 gcm~ 3 . 

Since the central star temperature drops down to less than 10 9 K within a 
few minutes after birth ||, the effects of finite temperatures on the equation 
of state can be neglected to a very good approximation. Consequently, the 
four differential equations (0)-(||) do not depend on time and thus need to be 
solved only once, at the beginning of the numerical cooling simulation. 

The thermal evolution is described by the partial differential equations of 
energy balance ||, 



d(Le^)_ 1 _ d(Te 



ON n\ dt 

and thermal energy transport, 



e v e^+Cy-^r- ^ , (13) 
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Table 1 

Input quantities for the cooling simulations performed in this paper 
Parameter References 

equation of state: 
crust 
core 



lay 

see table 



Superfluidity 



see table || 



Heat capacity 



I HP! 



Thermal conductivity: 
crust 
core 



\1Mb 



[11 18 1 



Neutrino Emissivity in the crust: 
pair-, photon-, plasma-processes 
bremsstrahlung [ 20 , 21 1 

in the core: see table || 

photosphere (no magnetic field) M 



d{Te 4 
dN 



16ir 2 r 4 Kn 



(14) 



The microphysical input into Equations ([13]) and ([H]) are the neutrino emis- 
sivity per unit volume, e„(P, T), heat capacity per unit volume, c v (P, T), and 
thermal conductivity, k(P, T). The boundary conditions for ( jl3|) and ( |I~4|) read 



L(iV = 0) = 0,, (15) 
T(N = N m ) = T m (r mj L m} M m ) } (16) 

where T m (r m , L m , M m ) is fixed by the properties of the photosphere at r = r m 
|7]|| . The initial condition for the star's temperature can be assumed, without 
loss of generality, to be T(r) = 10 11 K, because the star's thermal evolution 
at times greater than say a few months does not depend on the exact initial 
temperature profile inside the star. Equations ( p~3| ) and flU]) were solved nu- 
merically by means of a Newton- Raphson-like algorithm. Table p] summarizes 
the microphysical input quantities which enter the set of thermal structure 
equations, (|13|) and (|14D. The input quantities will be discussed in 

detail in sections 
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3 Models for the equation of state 



3.1 General remarks 



The equation of state (i.e., pressure and composition as function of energy den- 
sity) of neutron star matter is the basic input quantity whose knowledge over 
a wide range of densities, ranging from the density of iron at the star's surface 
up to about 15 times the density of normal nuclear matter reached in the cores 
of the most massive star of a sequence, is necessary when solving the ther- 
mal structure equations. Since one is dealing with highly isospin-asymmetric, 
net-strangeness carrying matter in /3-equilibrium whose properties cannot be 
explored in laboratory experiments, one is left with models for the equation of 
state which depend on theoretically motivated assumptions and/or specula- 
tions about the behavior of superdense matter. Sources of uncertainty concern, 
for instance, a competition between non-relativistic Schroedinger-based treat- 
ments versus relativistic field theoretical ones; moreover there exists even a 
considerable uncertainty with respect to the adequate many-body technique, 
which is to be introduced when solving the coupled equations of motion of 
many-baryon matter. Other items of uncertainty concern the baryon/meson 
composition of neutron star matter. This question has been treated in the liter- 
ature at different levels of complexity. The simplest description approximates 
neutron star matter by pure neutron matter, which however is certainly not 
the ground state of neutron star matter. In fact neutron matter will quickly 
decay by means of the weak force into chemically equilibrated neutron star 
matter, whose fundamental constituents-besides neutrons-are protons, hy- 
perons and possibly more massive baryons that become populated up to the 
highest densities reached in the cores of neutron stars. Finally we mention the 
possible transition of confined baryonic matter into quark matter, which, as 
demonstrated recently by one of us |22| , may set in at a density less than twice 
normal nuclear matter density. Such densities are easily reached in neutron 
stars possessing canonical masses of 1.4M . Last but not least, there is the 
open problem concerning K- or 7r-meson condensates in neutron stars. 



The cross section of a neutron star can be split roughly into three distinct 
regimes. The first one is the star's outer crust, which consists of a lattice of 
atomic nuclei and a Fermi liquid of relativistic, degenerate electrons. The outer 
crust envelopes what is called the inner crust, which extends from neutron 
drip density, p = 4.3 x 10 11 gcm~ 3 , to a transition density of about p tr = 
1.7 x 10 14 g cm -3 ||23|| . Beyond p tT one enters the star's third regime, that is, its 
core where all atomic nuclei have dissolved into their constituents, protons and 
neutrons. Furthermore, as outlined just above, due to the high Fermi pressure 
the core will contain hyperons, eventually more massive baryon resonances, 
and possibly a gas of free up, down and strange quarks. Finally ir- and K- 
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meson condensates may be found there too. The equation of state of the outer 
and inner crust has been studied in several investigations and is rather well 
known. We shall adopt the models derived by Baym, Pethick and Sutherland 
10 1 and Negele and Vautherin [11], respectively, for it. The models for the 



equation of state of the star's core will be discussed in detail in the next 
section. They fall into two categories, non-relativistic variational equations of 
state and the relativistic field theoretical ones. 



3.2 Non-relativistic variational approximation 



The non-relativistic models for the equation of state that are considered in our 
investigation use modern two-body interactions, Vy , whose parameters are ad- 
justed to the properties of the deuteron and the free two-nucleon scattering 
problem. These two-body interactions are supplemented with phenomenologi- 
cal three-body interactions, V^. The additional free parameters are adjusted 
to the properties of nuclear matter and the properties of 3 H and 4 He . The 
Hamiltonian's structure is given by 



H = E( V 4 2 + £ Vy + £ V ijk . (17) 

' i<j i<j<k 



The many-body method adopted to solve the many-body Schroedinger equa- 
tion is the variational approach [^4j], where a variational trial wave-function, 
\ty v >, is constructed from a symmetrized product of two-body correlation 
operators, F^, which act on an unperturbed ground-state, i.e., 



n f ^ 



|$ > . (18) 



The quantity |$ > denotes the antisymmetrized Fermi-gas wave function. 
The correlation operator contains variational parameters which are varied to 



minimize the energy per baryon for a given baryon number density n [24 





H 




< ^ 


V 


$v > 



E v (n) = min , ' " . (19) 



Equations fllBD and ([19]) were solved in P^] simultaneously in combination 
with the equations for electric charge neutrality and /^-equilibrium (q B and 
fi B denote electric charge and chemical potential of baryon B, respectively; 
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/i A is the chemical potential of leptons, A = e , /i ), 

E + l)|f - £ |f = , (20) 

B=p,n,... L=e,fj, 

fi B = fi n -q Bf i e , /i^ = /i e , (21) 
which relate the Fermi momenta, fcpB and &f,a and chemical potentials of the 



various baryons and leptons with each other [ 2l| . 

As representative models for the equation of state derived within this frame- 



work, we have selected a model derived by Wiringa, Fiks and Fabrocini |24 
where the two-body interaction UVi 4 supplemented with the three-body in- 
teraction UVII, derived by the Urbana group, serve as an input. This equation 
of state takes only protons and neutrons in equilibrium with respect to the 
weak force into account. 



3.2.1 Relativistic, field theoretical approximations 

The structure of a lagrangian that accounts for the physics of neutron star 
matter is given by |25| 



£ = E £°b + E £°m + E £°a + E 4 n L + , (22) 



^int 

B M A B,M 



where B is summed over those baryons that become populated up to the 
highest densities reached in the cores of neutron stars, i.e., B = n, p, E^ , A, 
H 0, ~, A~. The quantity M sums the mesons that are exchanged between them, 
M = a, u>, ii, p, rj, 5, and 0. The only contributing leptons are electrons and 
muons, hence A = e~, ur . The first three terms in (p2[) are the lagrangians of 
free particles. The term C l ^\j describes the interactions between the baryons, 
and Z^ " 4 -* represents self-interactions of a mesons. 

By means of applying the Green's function technique to the equations of mo- 
tion that follow from (^), one obtains the following coupled system of equa- 



tions p5|-[28[]. The first one is Dyson's equation which determines the two-point 
Green's functions, gf , 

(G#V(1, 2) - S B (1, 2)) gf(l, 1') = 5(1, 1') , (23) 



The second equation determines the effective baryon-baryon scattering ampli- 
tude in matter, T BB , 
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(12\T bb '\1'2') = (12\v bb '\1'2' -2'1') 



BB'\llnl o'i'\ 



+ (12|^|34)A^ (34,56)(56|T^^ |1'2') 



(24) 



where A stands for the Brueckner propagator [25,27-30 



The system of equations is closed by the expression for the self energy of a 
baryon in matter, S B , 



S B (1,2) 



t(13\T BB 142)0! (4, 3) . 



(25) 



We use the convention to sum or integrate over doubly occurring variables. 
The one-boson-exchange potential, v, sums the contributions arising from the 
various types of mesons, 

(12\v BB '\34) = £ (12|< B '|34). (26) 

M=a,uj,ir,p,... 



To obtain a relativistic model for the equation of state of neutron star matter, 
equations (|23| ) to are to be solved self-consistently in combination with 
relations (]20f) and ([H]) for electric charge neutrality and chemical equilibrium. 
This has been done for different levels of approximation. With increasing level 
of complexity, these are the relativistic Hartree (RH), relativistic Hartree- 
Fock (RHF), and relativistic Brueckner-Hartree-Fock (RBHF) approximations 
|25| , |31"lj . The latter scheme constitutes the most complex treatment that is 
presently manageable. It has been applied successfully to the composition of 
neutron star matter only recently. In contrast to RH and RHF, where the 
parameters of the theory, are to be adjusted to the properties of infinite nu- 
clear matter, the RBHF method makes use of one-boson-exchange potentials 
whose parameters are adjusted to the properties of the deuteron and the rel- 
ativistic nucleon-nucleon scattering problem. This leaves one with a so-called 
"parameter-free" treatment. 



3. 3 Comparison of the equations of state 



An overview of the collection of equations of state used in this paper is given 
in table 0. This collection serves to study both the influence of competing 
many-body approximations as well as different models for the nuclear forces 
on the cooling behavior of neutron stars. Table [| contains the corresponding 
nuclear matter properties. A graphical illustration of some of the equations of 
state is given in fig. [I]. 

With the exception of RBHF(B), the coupling constants of all other relativistic 
approximation schemes are chosen such that the properties of infinite nuclear 
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Table 2 

Dynamics and approximation schemes for equations of state derived for the cores 
of neutron stars 



EOS 


Composition 


Interaction 


Many body ap- 
proach 


Reference 


TJV-i/i+UVTT 




2- nuclci potential Urbana V14 and 

3- nuclei potential Urbana VII 


non-relativistic 

variational 

method 


l|f_3 


H V 


n n A V±>° CO,- „- 


Tl Vl d T1 (TO /"\T /T 'it /I ItlOtVK) c 

UjXdletlltiL, Ol L* , UJ , U IliesuilS 


RH 




RBHF(B)+HFV 


p, n, A, S*' , E°~, 
A~, e~, 


Exchange of cr, cxj, 7t, 77, (5 mesons 


RBHF 


¥ ' IF Ml 


G300 


p,n,A, Ei' ,^ '-^-, 


Exchange of cr, ui, p mesons 


RH 




^300 


p,n,A, E^ ,- '-^", 
/i — , pion-condensate 


Exchange of cr, ui, p mesons 


RH 




r^K240 
^13180 


p, n, A, S±.°, 3°.", 
e — , fi~ , quark matter 
phase 


Exchange of cr, u), p mesons 


RH, MIT bag 
model 




^K240 


p,n,A, S ± ' ,E°-,e-, 


Exchange of cr, ui, p mesons 


RH 





Abbreviations: RH=relativistic Hartree approximation, RBHF=relativistic 
Brueckner-Hartree-Fock approximation. 



Table 3 

Properties of nuclear matter for the equations of state used in this work 



EOS E/A [MeV] n [far 3 ] K [MeV] M* a SJ [MeV] 



non-relativistic mod- 
els 












UVu+UVII 


-11.5 


0.175 


202 


0.79 


29.3 



relativistic models 



HV 


-15.98 


0.145 


285 


0.77 


36.8 


RBHF(B) 


-15.73 


0.172 


249 


0.73 


34.3 


G300 


-16.3 


0.153 


300 


0.78 


32.5 


^300 


-16.3 


0.153 


300 


0.78 


32.5 


pK240 
^B180 


-16.3 


0.153 


240 


0.78 


32.5 


pK240 
^"M78 


-16.3 


0.153 


240 


0.78 


32.5 



The entries are: saturation density, no; binding energy per baryon, E/A; 
compression modulus, K; effective mass in units of the nucleon mass, M* = 
m* /m; asymmetry energy, agy. 
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Fig. 1. Pressure-density relation for several EOSs. The dots mark the highest central 
mass density reached in the most massive neutron star of each sequence. Neutron 
stars with still higher central densities are unstable against compressional modes. 

matter (n « (0.15-0.16) fm~ 3 , E/A « -16 MeV, m* w 0.78, K = 200-300 
MeV, a sy « 32 MeV) can be reproduced. For RBHF(B), Brockmann's one- 
boson-exchange potential (version 'B') was used as an input, which, as pointed 
out above, leads to a so-called parameter-free treatment. 

The ratio of hyperon to nucleon couplings to the meson fields, 

Xa = 9Ha/9Na, %u = 9Kw / 9Nlu , X p = g H p/gNp, (27) 



are not defined by the ground-state properties of normal nuclear matter and 



so must be chosen according to other considerations ||35|| . The binding energy 
of the A, for example, serves to strictly correlate the values of x a and x u , 
but leaves a continuous pairwise ambiguity which hypernuclear levels and 
neutron-star masses are able to resolve In this paper we choose for the 



hyperon-to- nucleon scalar coupling x a a minimum allowable value of x a ~ 0.5 
found from the lower bound on the maximum neutron-star mass that is also 
consistent with the A binding in nuclear matter. The corresponding value of 
x u amounts ~ 0.52. 

The possible transition of confined hadronic matter into quark matter at high 
densities is contained in the equation of state labeled Gj^fo J34|. The transi- 
tion was determined for a bag constant of B 1 / 4 = 180 MeV which places the 
energy per baryon of strange quark matter at 1100 MeV, well above the energy 
per nucleon in infinite nuclear matter as well as the most stable nucleus, 56 Fe 
(E/A ~ 930 MeV). (The possibility of absolutely stable strange quark matter 
will be considered in sect. f§.) Most interestingly, this model predicts a transi- 
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10 1 
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§10 13 



10 



10' 



10' 




r [km] 

Fig. 2. Density profiles of neutron star with mass 1.0 M constructed for a sample 
of different EOSs. 



tion to quark matter at a density as low as 1.6 n . The pure quark-matter phase 
is reached at a density of about 10 no, which is close to the central density of 
the heaviest and thus most compact star constructed from such an equation of 
state. We stress that these density thresholds are rather different from those 
obtained by numerous other authors, who performed work on this topic prior 
to Glendenning's paper. The reason for this lies in the seemingly innocuous 
assumption of treating neutron star matter as a simple body, that is, as matter 
that is characterized by only one conserved charge. This is not true for neu- 
tron star matter which possesses two conserved charges, baryonic as well as 
electric. Correspondingly, there are two chemical potentials rather than one, 
and the transition of baryon matter to quark matter is to be determined in 
the three-dimensional space spanned by pressure and the chemical potentials 
of the electrons and neutrons. The only existing investigations which accounts 
for this properly are those of Glendenning [E^,R3,RH] . 



The stiffness of the equation of state depends strongly on the internal degrees 
of freedom. Generally, the more degrees of freedom accessible to a system, 
the softer the equation of state. A softer equation of state, in turn, leads to 
neutron stars that are less massive but more dense than those obtained for 
stiffer equations of state. This is illustrated in figs. |^ to ^, where the density 
profiles of neutron stars constructed for different equations of state are shown. 
The stiffness of these equations of state can be inferred from fig. |I[ 
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Fig. 3. Same as fig. || but for a neutron star with mass 1.4 M & . 
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Fig. 4. Same as fig. [2] but for the heaviest neutron star of each sequence. 
4 Strange Stars 



The hypothesis that strange quark matter may be the absolute ground state 
of the strong interaction (not 56 Fe) has been raised independently by Bodmer 
57 1 and Witten 3£| . If the hypothesis is true, then separate classes of compact 



stars could exist, which range from dense strange stars to strange dwarfs to 
strange planets p9|-4Tl[] . They form distinct and disconnected branches of com- 
pact stars and are not part of the continuum of equilibrium configurations that 
include ordinary white dwarfs and neutron stars. In principle both strange and 
neutron stars could exist. However if strange stars exist, the galaxy is likely 
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to be contaminated by strange quark nuggets which would convert all neutron 



stars that they come into contact with to strange stars |42|-p4l. This would 
imply that the objects known to astronomers as pulsars are probably rotat- 
ing strange matter stars and not neutron matter stars, as is usually assumed. 
Unfortunately, lattice QCD calculations will not be accurate enough in the 
foreseeable future to give a definitive prediction as to the absolute stability of 
strange matter so that even to the present day there is not sound scientific ba- 
sis on which one can either confirm or reject the hypothesis, so that it remains 
a serious possibility of fundamental significance for compact stars plus various 



other physical phenomena |45j . Below we will explore the cooling behavior of 
strange stars and compare it with the one of their non-strange counterparts, 
i.e., ordinary neutron stars. This enables one to test the possible existence of 
strange stars and thus draw definitive conclusions about the true ground-state 
of strongly interacting matter. 



As pointed out by Alcock, Farhi, and Olinto |j46| , a strange star can carry a 



solid nuclear crust whose density at its base is strictly limited by neutron drip. 
This is made possible by the displacement of electron at the surface of strange 
matter, which leads to a strong electric dipole layer there. It is sufficiently 
strong to hold open a gap between ordinary atomic (crust) matter and the 
quark-matter surface, which prevents a conversion of ordinary atomic mat- 
ter into the assumed lower-lying ground state of strange matter. Obviously, 
free neutrons, being electrically charge neutral, cannot exist in the crust, be- 
cause they do not feel the Coulomb barrier and thus would gravitate toward 
the strange-quark matter core, where they are converted by hypothesis into 
strange matter. Consequently, the density at the base of the crust (inner crust 
density) will always be smaller than neutron drip density. The situation is 
illustrated in fig. |5| where the equation of state of a strange star with crust is 
shown. 

The strange-star models presented in the following are constructed for an 



equation of state of strange matter derived by Farhi and Jaffe | 47| . The value 
chosen for the bag constant and strange quark mass are B 1 / 4 = 145 MeV and 
m s = 150 MeV, respectively. These choices place the energy per baryon of 
strange quark matter at 850 MeV, which corresponds to strongly self-bound 
strange matter. 



5 Superfluidity 



The effective nucleon-nucleon interaction contains a rather strong repulsive, 
short-range part and a weaker, attractive long-range part. Hence, depending 
on the average distance of the particles, nuclear matter may become super- 
fluid (see for a recent review). Neutrons, for instance, are believed to 
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P [9 cm" 3 ] 

Fig. 5. Pressure-density relation for strange stars with crust. The bag constant is 
B l / A = 145 MeV, the mass of the strange quark is m s = 150 MeV. The solid dot 
marks the central density reached in the heaviest stable strange star. Stars with still 
higher central densities are unstable against compressional modes. 



form 1 So superfluid pairs in the density range between neutron drip and and 
the saturation density of normal nuclear matter (see below). At supernuclear 
densities, i.e., in the core region of neutron stars, the repulsive interaction be- 
comes dominant and the 1 So superfluidity gap closes. However, the tensor and 
spin-orbit interaction become attractive at such densities and the formation of 
3 P2-neutron pairs becomes possible. Finally, the proton concentration in the 
core region may reach a partial density simlar to that of the pair-condensed 
neutron gas at subnuclear densities, enabling the protons to form superfluid 
1 S pairs. Since the protons are charged, the proton superfluid is a supercon- 



ductor. The different superfluid regimes can be summarized as follows |Ej -|51] : 



(i) 7 x 10 n gcm- 3 to 2 x 10 14 gcm" 

(ii) 2 x 10 14 gcm- 3 to 4 x 10 14 gcnr 
(hi) 2 x 10 14 gcm~ 3 to 5 x 10 16 gcm" 



neutrons in 1 S -pair state, 
protons in 1 S -pair state, 
neutrons in 3 P 2 -pair state. 



One has to keep in mind that these values depend on the composition of 
matter. The above values are computed for the equation of state labeled HV. A 
further possibility not mentioned so far are 3 D 2 neutron-proton pairs. Finally, 
it is possible that hyperons and quarks form superfluid states too |]52| , |53|| . Until 
now, however, no detailed calculations on this topic exist and one is thus left 



with parameter studies (cf. section |?\7|) . 



If protons and/or neutrons become superfluid, the neutrino emissivity, thermal 
conductivity, and heat capacity are reduced by an exponential factor of the 
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Table 4 

Gap energies and critical temperatures of superfluid states in neutron star matter 



Neutron 1 Sr 



Neutron 3 P 2 Proton 1 Sq 



Reference 



m 







A max 



[MeV] 
[10 9 K] 



1.13 
7.4 



0.62 
0.8 



0.25 
1.6 



from p [gem 3 ] 7 x 10 11 
top [gem" 3 ] 2 x 10 14 



2 x 10 14 2 x 10 14 
5 x 10 16 4 x 10 14 



The density ranges are calculated for the HV equation of state. A^ ax = 
1.13 MeV is obtained for the 4 So superfluid gap as estimated by Ainsworth 
et al. Using the gaps estimated in the more recent paper by Wambach et 
al. leads to A^ ax = 1.0 MeV, which influences the cooling of neutron 
stars only insignificantly. 




10 

p [g cm" 3 ] 



Fig. 6. Energy gap versus mass density computed for the HV EOS. 



form[|]exp(— A sf /fc B T), where A sf is the gap energy |r7||56|. Unfortunately, the 
gap energies are not very well known because of uncertainties in the effective 
baryon masses and polarization effects []57|]. We will take the gap energies 



for protons from 
from ^{| and 



and those for the So and P 2 states of neutrons 
respectively. An overview of the gap energies and critical 
temperatures is put together in table As an example, we show in fig. § 
the density dependence of A sf for one of our model equations of state. A 
qualitatively similar behavior is obtained for the other models. 



For the 3 P2 state this is rigorously true only for isotropic matter |53|-55]. 
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Table 5 

Neutrino processes studied in this work, their associated neutrino emissivities, e u , 
and the critical threshold densities, n CT [ t , below which these processes are forbidden 

e u [erg s" 1 ] n crit [fm~ 3 ] 



Process 



Reference 



Nucleon bremsstrahlung 
modified Urea 
direct Urea 
Pion-condensate 
Kaon-condensate 
Quark /3-decay 



10 19 x T 9 8 
10 21 x T 9 8 
10 27 x Tl 
10 26 x Tl 
10 25 x Tl 
10 24 x Tl 



5§|58 



0.25 - 0.5 
0.21 
0.71 
0.29 



1 



The first two processes are slow ones. They are possible at all densities, 
for which bulk nuclear matter exisits. The critical density for the direct 
Urea process depends on the equation of state, which leads to a range of 
threshold densities. 



- Bremsstrahlung 




I modified Urea 




r direct Urea 




[ pion Urea 




kaon Urea 
























/ 


1 /' 




/ 

/ 

1 


1/' 1 









P [g/cm 3 ] 

Fig. 7. Neutrino emissivities of different neutrino-emitting processes as a function of 
energy density for T = 10 9 K. The particle composition is taken from the Gfgn EOS. 

6 Neutrino emissivities 



The neutrino emission processes can be divided into slow and fast ones. The 
main feature, which distinguishes fast from slow processes, is the existence of 
a certain threshold density for the fast ones, below which they cannot take 
place. An overview of all processes which we shall take into account is given 
in table |5[ The associated neutrino emissivities are shown in fig. 0. 
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6. 1 Neutrino pair Bremsstrahlung 



The most important neutrino process in the crust is neutrino-antineutrino pair 
production by electron bremsstrahlung. In early estimates it was found that 
the contribution of this process to cooling is competitive with standard cooling 
via the modified Urea process. However, recent investigations have shown that 
the mass of the ionic crust is reduced by a factor of two to three if one considers 
the recently established value p tr ~ 1.7 x 10 14 gcm -3 for the transition density 
between crust and core [23], instead of the older value p tT ~ 2.8 x 10 14 gcm -3 



||10|| . Additionally, considerations by Pethick and Thorsson [21J have pointed 



out that the standard perturbative treatment of determining e v is not suf- 
ficiently accurate below a temperature of 10 10 K. Furthermore taking the 
electronic band structure into account can reduce the energy emission rates 
exponentially for temperatures below 10 10 K. We include these corrections to 
the emissivities derived by Itoh et al. in p^,E5[. Additionally to the above 



process, pair, photon and plasmon processes may also be important in the low 



density regime. Their effects have been estimated too by Itoh et al. in |19 
and will be used in our calculations. 

Neutrino-antineutrino pair bremsstrahlung production can also occur in the 
core according to 



n + n— >n + n + i/ + P, (28) 
n + p^n + p + z/ + z/, (29) 
e~+p^e~+p + z/ + z/. (30) 



We note that the emissivity of the last bremsstrahlung process, Equation (p0|), 
becomes very important if the neutrons become superfluid. In this case this 



process dominates all other processes. For the processes (|28|) to (|30|) we use 
the emissivities as calculated by Friman and Maxwell [jSG^jSfifl. 



6.2 Modified Urea process 

The beta decay and electron capture processes, 

n — ► p + e~ + u e and p + e~ — > n + v e , (31) 



also known as the direct Urea processes, were not included in cooling calcu- 
lations p6H7(J until recently, since energy and momentum conservation can 



only be fulfilled if the proton fraction exceeds a critical value (cf. section |6.3| ) . 
Most interestingly, many of the older, non-relativistic models for the equations 
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of state of neutron star matter indeed do not lead to sufficiently high enough 
proton fractions. Consequently, the so-called modified Urea process, which is 
characterized by the occurrence of a bystander particle which is necessary to 
conserve both energy and momentum of the baryons that are being scattered, 
i.e., 

n + n— ► n + p + e~ + z? e and n + p + e~ n + n + z/ e . (32) 



were considered to be the dominant cooling mechanism of neutron stars. The 
emission rates associated with processes were calculated by Friman and 
Maxwell ||58|| . There the exchanged 7r meson is treated as a free particle. Such 
an approximation neglects the influence of the medium on the propagating 
pion. Furthermore, the exchanged pions undergo weak transitions (so-called 
exchange current effect). An estimate of the influence of such effects, performed 
by Voskresenskii and Senatorov [ffll , [72| , shows that the neutrino emissivities 



may increase by up to several orders of magnitude. The implications of these 
effects of cooling will be studied in a different investigation 173 . 



6.3 Direct Urea process 



As already mentioned, the direct Urea process (|31| ) is only possible if the 



proton fraction exceeds a certain critical value |59],[71[. If the star is assumed 



to be made up of only neutrons, protons and electrons, the critical proton 
fraction amounts about 11%. In one includes hyperons and muons too, the 
critical value increases slightly to 13%. In contrast to older, non-relativistic 
equations of state which predict proton fractions that lie below the critical 
values, modern non-relativistic as well as relativistic equations of state contain 
sufficiently large enough proton fractions such that the direct Urea process 
becomes possible. The proton fraction depends crucially on the asymmetry 
energy which unfortunately is not known at higher densities. From our study 
we find that the proton fraction is generally larger for the relativistic equations 
of state than for the non-relativistic ones. This trend can be seen in fig. [8], where 
the proton fractions are shown for several equations of state of our collection. 
The critical energy densities beyond which the direct Urea becomes possible, 
are indicated. Because of the different asymmetry behaviors of these equations 
of state, the critical densities vary considerably from one equation of state 
to another. This, in turn, links the onset of the direct Urea process to the 
condition of a sufficiently large enough star mass. For the emissivities of the 
direct Urea process we shall use those calculated by Lattimer et al. J59|. The 



emissivity of the direct process exceeds the one of the modified process by up 
to ten orders of magnitude, at low temperatures of T « 10 7 K. 

So far we have considered the direct Urea process in neutron star matter made 
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10 




10' 



Fig. 8. Proton fraction in neutron star matter as a function of energy density for 
different equations of state. The approximate critical proton fraction n p /nB = 0.13 
for the direct Urea process is marked by the horizontal line. The exact critical energy 
densities together with the corresponding star (masses in units of M©) are marked 
by arrows. 

up of protons and neutrons only. Strange baryons and A isobars, however, 
which may exist in such matter rather abundantly too, also produce neutrinos 
via direct Urea processes of the following type, 



A + e + v e and A + e — > £ + v c 



(33) 



Such processes have been studied first by Prakash et al. |75| , and the implica- 
tions for cooling of neutron stars were explored in [ 76| . We have included the 
above processes in cooling simulations elsewhere |[77|| . It was found that the 
cooling behavior due to direct Urea among the nucleons only dominates over 
the energy loss associated with processes the (|3"3]). Or, in quantitative terms, 
the cooling curves obtained by restriction to the direct Urea process among 
non-superfluid nucleons (section |7.4| ) remain practically unchanged. 



6.4 Quasi-particle processes in boson-condensed matter 



In the charged pion (n c -) condensed phase, the quasi-particle Urea processes 



77(71") — > 77(71") + e + z/ e and 77(71") + e — > 77(71-) + v e (34) 
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are possible. The quasi-particle state r\ is given by the superposition of proton 
and neutron states, 



\v±) = cos 6>|n) =1= i sin #|p) 



(35) 



with the mixing angle 9. The emissivity for this process was calculated by 
Maxwell et al. [0. 

In recent years, some authors made a strong case for the condensation of kaons 
in the cores of neutron stars (see, for instance, |78|-p0|). The critical density 
beyond which this process may take place is uncertain. Estimates lie in the 
range of 3 — 5n [|3l] , |31] ,|3^|. The emissivities of the quasi-particle direct Urea 
processes, 

n(K) -> p(K) + e" + v e and p(K) + e - -> n(K) + u e , (36) 



are reduced, relative to the uncondensed case, by a factor of cos 2 (9/2) fl62|j . 
Additionally to (|36|) there are two more processes, 



n(K) — > n(K) + e + u e and 
p(K) — > p(K) + e~ + u e and 



n(K) + e~ -> n(K) + u e , 
p(K) + e~ -> p(K) + u e , 



(37) 
(38) 



whose emissivities were calculated by Thorsson et al. [62]. The density depen- 
dence of the mixing angle 9 and the value of the critical density, n c = 4.18 n 



are taken from 61 . 



6. 5 Quark processes 



The possible existence of quark matter in the cores of neutron stars has already 



been outlined in sect. O. The latest calculations on this topic predict an onset 
of the transition of confined hadronic matter into quark matter at n Pd 1.6 no- 
(This is not to be confused with the possibility of absolutely stable strange 
matter, to which we will return in sect. HI). 



The following weak processes are important in quark matter consisting of up, 
down and strange quarks, 



d — > u + e + v e and u + e 
s — > u + e~ + u e and u + e~ 



d + v e , 
s + v c . 



(39) 
(40) 



We shall use the neutrino emission rates given by Duncan, Shapiro, and 
Wasserman pBJ to simulate the cooling history of neutron stars containing a 



21 



— our calculation 
■ Van Riper (1991) 




Fig. 9. Comparison of the computed surface temperature versus age with the results 



of Van Riper [12|. The parameters are: M = 1.4M©, BPS EOS [1G|, no superfluidity 
and no enhanced cooling processes are included. 

mixed phase of baryons and quarks in their interior (hybrid stars) and strange 
stars. 



7 Results and comparison with observed data 



7. 1 Test of the numerical code 



We tested the correctness of our numerical code by means of comparing the 
outcome of several cooling calculations with those obtained by other authors. 
As an example, in fig. |9] we compare our surface temperatures, obtained for 
a particular model for the equation of state, with the corresponding results 
obtained by Van Riper [|T!J. The agreement is extremely good for all star 
ages. To achieve this agreement, however, we had to reduce the core's thermal 
conductivity, which was not specified by Van Riper, from the value given in 
fP7| by a factor of 40. Otherwise deviations of about ~ 20% were obtained at 
star ages between r ~ 1 to 10 years. The perfect agreement at later stages 
however is not influenced by this. 
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Table 6 

Survey of observed pulsar ages and luminosities, denoted r and L, respectively 



Pulsar 



Name 



logr [yr] logL [erg/s] Reference 



not enough data available for spectral analysis 
4.25 32.8 ±0.7 
4.50 33.2 ±0.6 
4.61 33.1 ±0.4 



1706-44 
1823-13 
2334±61 



83 
84 
M 



power-law-type spectra or spectra with only a high temperature component 



0531±21 


Crab 


3.09 


35.5 ±0.3 


1509-58 


SNR MSH 15-52 


3.19 


33.6 ±0.4 


0540-69 




3.22 


36.2 ±0.2 


1951+32 


SNR CTB 80 


5.02 


33.8 ±0.5 


1929±10 




6.49 


28.9 ±0.5 


0950±08 




7.24 


29.6 ± 1.0 


J0437-47 




8.88 


30.6 ±0.4 



86 



89 



urn 



92 



m 



spectrum dominated by a soft component 
0833-45 Vela 4.05 32.9 ± 0.2 

0656±14 5.04 32.6 ± 0.3 

0630±18 Geminga 5.51 31.8 ± 0.4 

1055-52 5.73 33.0 ± 0.6 



Ml 



95 
9§ 
071 



The age is determined by the spin down of the pulsar, the luminosities by 
spectral fits or the totally detected photon flux, depending on the categories 
(i) to (iii) listed in the text. 



7.2 Body of observed cooling data 



From the X-ray observatories Einstein, EXOSAT, and ROSAT the X-ray spec- 
tra of at least 14 pulsars have been observed to date. Unfortunately, the data 
are not always sufficient to extract the effective surface temperature of the 
corresponding pulsar. One can classify the spectra into three categories |4j] 
(compare with table |6|): 

(i) The photon yields of three pulsars (PSR's 1706-44, 1823-13, and 2334+61) 
contain too few photons in order to perform spectral fits. The luminosities 
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are calculated by using the totally detected photon flux. These pulsars 
are marked with triangles in the figures. 

(ii) The spectra of seven pulsars, including the crab pulsar (PSR 0531+21), 
can only be fitted by a power-law type spectrum or a blackbody spectrum 
with a very high effective temperature, and an effective area much smaller 
than the neutron star surface. Their X-ray emission is predominated by 
magnetospheric emission. Therefore, their temperatures, determined from 
the spectral fits, are probably too high. Theses pulsars are marked with 
solid dots. 

(iii) Finally, there are four pulsars (i.e., 0833-45 (Vela), 0656+14, 0630+18 
(Ge-minga), and 1055-52) that allow two-component spectral fits, where 
the softer blackbody components corresponds to the actual surface emis- 
sion of the neutron star and the harder to some magnetospheric emission. 
These pulsars are marked with squares. 

When discussing spectra of neutron stars one should keep in mind that neutron 
stars are not black bodies, because hydrogen and helium in their atmospheres 
modify the blackbody spectrum. Besides, strong magnetic fields can also af- 
fect the surface emission. With respect to the luminosity determinations we 
note that the luminosities determined via hydrogen- or helium-atmospheric 
models are generally lower by a factor of about three than those determined 
via blackbody spectral fits. The inverse is true for the temperature which is 
smaller by a factor of about two. An uncertainty in temperature by a factor of 
two varies the luminosity, according to the Stefan-Boltzmann law (L oc T 4 ), 
by a factor of 16. It is thus not appropriate to perform a comparison with sur- 
face temperatures of stars M. Instead one should be using the luminosities, 
as is the case here. 



7.3 Standard cooling calculations 



Figure |10| shows the cooling history of a neutron star with M = 1.4M that 
cools down via the standard scenario only, i.e., fast cooling processes are not 
allowed. The different curves correspond to different models for the equation 
of state. One sees that the standard cooling behavior depends only weakly 
on the equation of state. The UV14+UVII model cools somewhat faster at 
later times than the other two models. This is due to the rather soft behavior 
of the UV14+UVII equation of state at intermediate nuclear densities (fig. 
[I]), which leads to neutron stars with small radii (figs. @-@). Obviously, most 
of the observed data are not in agreement with the assumption of standard 
cooling. Specifically Vela's (0833-45) luminosity lies considerably below the 
computed curves. The data of PSR 0630+18 (Geminga) and a few others are 
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Fig. 10. Standard cooling behavior of a neutron star with mass M = 1.4M Q for 
several different models for the EOS. 



in relatively good agreement, however. 



7.4 Enhanced cooling via direct Urea process 



As mentioned above, one way by means of which neutron stars can cool very 



efficiently is the direct Urea process |98 H100| . This is demonstrated graphically 



in figs. [ll|HT3|for the equations of state UV14+UVII, HV, and Gsoo- The proton 
fraction in the stars with M = 1.0M Q is for all three cases below the threshold 
density for the direct Urea process. Hence these stars cool via the standard 
mechanism, the modified Urea process. The cooling behavior of all other stars 
at early stages is dominated by the direct Urea process. The sudden drop 
in luminosity at ages of several hundred years has its origin in the fact that 
the core cools down much faster through the enhanced cooling process than 
the crust. This causes a temperature inversion in young stars. After several 
hundred years, depending on the thickness of the crust, the 11 cooling wave" , 
built up by the temperature gradient, reaches the star's surface which leads to 
the sudden luminosity drops. (An estimate of the time needed by the cooling 



2 Note that there is a considerable uncertainty in the ages of pulsars that are derived 
from the braking index, n, which, due to lack of information, usually is taken to 
be constant during a pulsar's spin-down. Moreover a canonical value of n = 3 (T^j, 
which corresponds to a slow-down by emission of pure magnetic dipole radiation, is 
assumed. This too constitutes a rather crude approximation. 
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Fig. 11. Enhanced cooling of neutron stars via direct Urea process. The underlying 



EOS is UV14+UVII. The observed data are labeled in fig. 10 
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Fig. 12. Same as figure 11 but for HV. The observed data are labeled in fig. 10 



wave to reach the surface can be found in [|101|| .) The drop of luminosity 
increases with massive, since the region inside the core of a star where the 
direct Urea process is possible grows with density. The drop is most strongly 
pronounced for UV14+UVII. This has its origin in the 3 P2 neutron-superfluid 
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Fig. 13. Same as fig. 11 but for G300 The observed data are labeled in figure [[(]. 
which does not extent into the central core. 

Changes in the hyperon couplings strengths in matter influence the cooling 
curves only very little, as can be seen in fig. [14| for neutron stars of several 
different masses. (Note that this may be drastically different if the direct Urea 
process among hyperons occurs is dense stellar matter |75| , |76| .) The star with 
M = 0.7 M Q does not contain hyperons, since its core is not dense enough 
to reach the threshold density for hyperon population. Hence there is only 
one cooling curve in this case. The situation is different for the heavier stars 
shown in this figure, for which the hyperon concentration increases with star 
mass ||105|| . The hyperons are coupled to the meson fields either with the same 
strengths as the nucleons, i.e., x a = x^ = x p = 1 (cf. Eq. (p7|)), or are chosen 
such that mutual compatibility between the A binding energy in saturated 
nuclear matter, neutron-star masses, and experimental data on hypernuclear 
levels is achieved (cf. sect. |3.3|) . As for the other figures of this section, the 
direct Urea process among the nucleons is included, and no superfluidity is 
taken into account. 



7. 5 Enhanced cooling via pion condensation 



The impact of pion condensation on the cooling of neutron stars ||100|1 is 
demonstrated in fig. [TJ|. comparison with fig. [13| shows that the cooling be- 
havior of such stars is qualitatively similar to the one of stars which cool via 
the direct Urea process. The threshold density of the pion condensate lies at 
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Fig. 14. Enhanced cooling of neutron stars due to direct Urea among nucleons, but 
for varying hyperon coupling strengths (see text). The underlying EOS is Gj^g . 
The observed data are labeled in fig. HI]. 
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Fig. 15. Enhanced cooling of neutron stars due to pion condensation. The underlying 
EOS is GJqq. The observed data are labeled in fig. [[(]. 
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Fig. 16. Enhanced cooling of neutron stars due to kaon condensation. The underlying 
EOS is Gsoo- The observed data are labeled in fig. [u]. 

3.6 x 10 14 gem -3 , just slightly above normal nuclear matter saturation den- 



sity. Hence all stars displayed in fig. |D| with the exception of M = O.5M , 
possess a pion condensate and therefore show enhanced cooling triggered by 
the condensate. 



7.6 Enhanced cooling via kaon condensation 



The cooling behavior of kaon-condensed stars (see also p3| , |100| , ^02| - |104|| ) is 
qualitatively similar to those of pion-condensed stars. This is graphically il- 
lustrated in fig. pi] for the G300 equation of state. There is however a significant 
quantitative difference between both types of condensation which concerns the 
threshold densities. It amounts for kaon condensation about 4.8 n |61|, which 
is by a factor of ~ 2.5 larger than the threshold density established for pion 
condensation, ~ 1.3 no [|105|| . A density of 4.8 n is reached only in the most 
massive star model, M = 1.78 M , constructed for Gsoo- Hence the other two 
models cool via the modified Urea process only. For the purpose of compari- 
son, the central density of the M = 1.4 M & neutron star amounts 2.4 no which 
lies almost by a factor of two below the kaon threshold. The situation would 
be different for softer models for the equation of state, which lead to higher 
central neutron star densities. For the broad collection of modern equations of 



state of ref. [25| we find central densities in M ~ 1.4 M Q stars that range from 



about 2.3 n to 4.2 n . The latter value comes much closer to the critical value 
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Fig. 17. Enhanced cooling of neutron stars due to an admixture of quark matter 
in their cores (G^gg EOS). Curves corresponding to the G^||° EOS are without 
quarks. Enhanced cooling via direct Urea is taken into account too. The observed 
data are labeled in fig. Il0]. 



for kaon condensation derived in |61| , which, of course, is model dependent 



too and thus not accurately known. Clearly, the smaller the threshold density 
the faster the cooling of compact stars via kaon condensation will be. 



7.7 Enhanced cooling due to the presence of quark matter in the cores of 
neutron stars 



As discussed before, modern investigations of the possible transition of baryon 
matter into quark matter predict transition densities less than twice normal 
nuclear matter density. Such densities are easily reached in neutron stars with 
masses M « 1.4M (see also discussion at the end of section |7T6| ). To 



simulate the cooling history of such (hybrid) stars we apply Glendenning's 
Gmlo equation of state, whose specific properties have been discussed in con- 
nection with table || The weak reaction processes of the quarks are those 
given in Equations (|3^) and ([4(]). They speed up the energy loss of a star 
considerably, as is illustrated in fig. fL7|, where the solid curves correspond to 
enhanced cooling caused by the quarks (see also [0). The dashed curves re- 
fer to stars with no quarks in their cores, which therefore cool more slowly. 
Note that the density in the light stars with M < 1 M is not high enough to 
reach the transition into quark matter. Enhanced quark cooling is therefore 
not possible for these stars. Finally, if the quarks should form a superfmid, 
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Fig. 18. Cooling behavior of strange quark matter stars of masses M = 0.7 M & and 
M = 1.85 Mq with and without crust. The quarks are treated as non-superfiuid. 
The observed data are labeled in fig. [H]. 

the emission of neutrinos is suppressed and cooling proceeds at a much slower 
pace (dotted curves). Bailin and Love have stated that quarks might become 
superfluid ||52|| . Unfortunately, there exists - as far as we know - no estimation 
of the gap energy of superfluid quark matter. Therefore, in order to simulated 
this effect, we simply assumed a plausible value of A sf = 0.1 MeV for the 
1 So pairing state of quarks, which was taken to be constant over the whole 
density range. 



7.8 Strange quark matter stars 



We restrict ourself to determine the cooling behavior of strange stars which are 
either bare or possess the densest possible hadronic crust that can be carried 
by such objects (inner crust density equal to neutron drip density). Figure p 



reveals that, no matter how thick the crust, such stars cool very rapidly ||106 



Among the body of observed luminosities, there is presently one data point 
that would be consistent with the assumption that the underlying object is 
a strange star. This interpretation is not stringent of course, because some of 
the enhanced cooling mechanisms studied above might do as well. If one treats 
the quarks as superfluid, neutrino production ceases and therefore cooling is 
considerably slowed down. To simulate this effect qualitatively we assumed a 
density- independent energy gap of A s f = 0.1 MeV. The numerical outcome 
of such a simulation is exhibited in fig. O. One sees that in this case the 
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Fig. 19. Cooling behavior of strange quark matter stars with crust. The masses are 
M = 0.7 Mq and M = 1.85 Mq. The quarks are treated as superfluid particles. The 
heat capacity is reduced by an exponential factor in one model, while in the other 
model it is kept unchanged. The observed data are labeled in fig. [l(]. 

cooling behavior of strange stars is similar to the enhanced cooling behavior 
of neutron stars (compare figs. pT|- [l7|) . Hence, with the exception of PSR 0540- 
69 practically all observed pulsar luminosities can now be reproduced too. A 
quantitative clarification of the question as to whether or not quark matter 
indeed becomes superfluid, and if so, at which densities, attains particular 
interest from astrophysics too. There is also the possibility that the quark 
flavors become superfluid at different densities. Consequently the heat capacity 
would be somewhat smaller. To simulate this possibility we reduce the heat 
capacity by an (arbitrary) factor of three. The consequences for the cooling 



behavior are shown by the dotted curves in 19 



In closing this section, we mention that one of the main differences between 
the cooling behavior of neutron stars and strange quark matter stars concerns 
the time needed for the cooling wave to reach the surface. This time is much 
shorter for strange quark matter stars because their crusts have a thickness of 
only ~ 0.2 km as opposed to ~ 2 km for stars with M ~ 1.4M . 



8 Conclusions 



We have already pointed out in section [7]2] that the observed cooling data of 
pulsars marked with squares in our figures appear to be the most significant 
ones. Primarily these data should be reproduced by theoretical cooling calcu- 
lations. As we have seen above (section |7.3|) the luminosities of a few pulsars, 
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Fig. 20. Influence of several different enhanced neutrino-emission processes on the 
cooling of neutron stars with masses M = 1.4 M©, except for the kaon-condensed 
model where the mass is M = 1.78 M & (The 1.4 Mq models do not possess large 
enough central densities to overcome the threshold density for kaon condensation). 
The observed data are labeled in fig. 1C. 



like 1951+32 and 1055-52, appear to be compatible with the standard cooling 
scenarios. Depending on the model employed for the equation of state, this 
may also be the case for Geminga and pulsars 0656+14, 2334+61 (cf. fig. ITUp. 



Others, however, like the famous Vela pulsar (PSR 0833-45) and PSR's 1706- 
44, 1929+10, which have rather low observed luminosities, may require the 
introduction of so-called fast (or enhanced) cooling processes, which are char- 
acterized by higher neutrino emission rates by means of which stars cool more 
rapidly. A final decision on this issue - standard versus enhanced cooling |3| 
is presently hampered by the uncertainties in a number of quantities, including 
the modified Urea neutrino-emission rate, superfluid gaps, equation of state, 
and mass and age of the star. Nevertheless, subject to the enormous progress 
that is being made presently in exploring the behavior of superdense matter 
p|, one may feel confident that this will change in the foreseeable future. 



Figure [2(] compares the influence of different enhanced neutrino emission pro- 
cesses, i.e., direct Urea, pion and kanon condensation, and quark Urea, on 
the cooling behavior of a neutron star with mass M = 1.4 M & . The nucleons 
and quarks are treated as non-superflid particles. One sees that the inclusion 
of any of these enhanced cooling processes reduces the star temperature too 
quickly in order to get agreement with the observed data points. The only ex- 
ception is PSR 1929+10. The most significant data points (squares), however, 
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lie considerably above the enhanced cooling curves. 



The enhanced cooling scenarios can be slowed down if one assumes that the 
neutrons in the cores of neutron stars are superfluid, as was the case in figs. 



1 I |u]. The only problem left then concerns the horizontal "plateaus" in these 
figures at star ages between 10 2 and 10 5 years which tend to lie somewhat too 
high for stars with canonical mass, M = 1.4 M Q , to be consistent with the 
data. Note, however, that one gets agreement with some of the data points 
if the star mass is assumed to be different from the canonical value. Another 
possibility to get agreement consists in varying the gap energy, A s f, which, as 
mentioned in section [5], is not very accurately know. This is particularly the 
case for the high-density 3 P 2 superfluid. A reduction of this gap by a factor 
of two, for instance, moves the cooling curves up right into the region where 



the oberved data is concentrated, as can be seen in fig. |21 
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Fig. 21. Influence of changes in the superfluid 3 P2 gap on the cooling of neutron 
stars constructed for Gsoo- The enhanced cooling processes are direct Urea in pion 
(solid curves) and kaon condensed (dotted curves) matter. The observed data are 
labeled in fig. |HJ 



Of course there are other possibilities, besides reducing the gap energy, by 
means of which agreement with the observed data may be achieved. For exam- 
ple, the superfluid phase of neutrons might not reach the center of very massive 
neutron stars, as is the case for stars constructed for UV14+UVII (see fig. 



Tip . This is due to the smaller proton fraction and the resulting higher Fermi 



momenta of neutrons. Another possibility is an additional cooling process not 
suppressed by superfluidity, as, for example, the superfluid pair-breaking pro- 
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cess [f71 ,fT3 1 , or so-called internal heating of neutron stars [|107H1H • 



9 Summary 



We have applied a broad collection of modern, field theoretical models for 
the equation of state of superdense matter to the investigation of the cooling 
behavior of compact stars, which comprise neutron stars and their strange 
counterparts, the hypothetical strange stars. The models for the equation of 
state were derived elsewhere, using many-body approximation techniques at 
different levels of sophistication and complexity. Furthermore, uncertainties in 
the behavior of matter at extreme densities, like the phenomenon of meson 
condensation, superfluidity, hyperon degrees of freedom, transition to quark 
matter, and absolute stability of strange quark matter are incorporated in this 
collection. 

Besides standard cooling, which, subject to uncertainties in the neutrino emission- 
rate of the modified Urea process, superfluid gaps etc. (cf. section ||), gives 
agreement between the theoretical cooling curves and the observed data for 
some but not all pulsars (the obtained luminosities are considerably too high 
in some cases), we studied various enhanced cooling mechanisms too. These 
are connected to higher neutrino-emission rates coming from the direct Urea 
process, 7r- and f^-meson condensates, and up, down and strange quarks in 
the cores of neutron stars. It turns out that the enhanced processes altogether 
cool the star models too rapidly in order to get agreement with the data. 
So one is left with introducing processes which delay somewhat the cooling 
as driven by the enhanced processes. A plausible candidate is superfluidity, 
of which there occur probably two different types in the cores of neutron 
stars, namely and 3 P2- Due to the complicated theoretical analysis, the 
density dependences and widths of the associated superfluid gaps cannot be 
determined exactly, especially in the latter case. By using the more reliable 
1 So gaps we find that already small changes of the theoretical value for the 
3 P2 gap shift the cooling curves into the region of observed data, giving an 
overall good agreement. An accurate determination of the superfluid gaps in 
neutron star matter seems therefore of particular interest. 

Most interestingly, we find that-for a given star mass-not a single model for the 
equation of state of our collection is able to explain more than just a few data 
points! This makes one wonder if the observed data are possibly associated 
with stars of different masses. We recall that the canonical neutron star mass 
is about 1.4 M , and the observed data ||112|| seem to concentrate about this 



value. Nevertheless, the existence of pulsars with quite different masses, e.g., 
neutron stars in binary systems that accrete mass from their companions, 
seems equally plausible. If one assumes that the observed luminosities are 
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indeed coming from stars with different masses, those of Vela and Crab could 
amount ~ 1 M Q and ~ 1.4 M & , respectively, then one may have a reason for 
the different cooling behaviors. 



As mentioned repeatedly, the absolute stability of strange quark matter would 
have implications of greatest importance for various branches of physics [^5 



It is agreed on that compact stars form ideal environments where such mat- 
ter could be found. Adopting a simple model for the EOS of strange matter, 
we studied the cooling history of stars made up of such matter, too. If the 
emissivities and heat capacities of quark matter used here constitute reliable 
estimates, which, at present, is an unresolved issue, we find that pulsars like 
Geminga, Monogem, and PSR 1055-52 could be ruled out as strange pul- 
sar candidates. There is however one pulsar, PSR 1929+10, which could be 
interpreted as a strange star. Of course, the QCD related uncertainties in 
the properties of strange matter do not allow to draw any stringent conclu- 
sions yet. In order to provide a feeling for how sensitive the cooling history of 
strange stars depends on the transport properties of strange matter, we also 
performed cooling simulations where the emissivity and heat capacity were 
slightly varied, due to a possible superfluid behavior of the quarks. As it turns 
out, this moves the cooling curves right into the region of the luminosity-age 
plane where most of the observed data are concentrated. So there is clearly a 
big need for reliably determined transport properties of strange quark matter. 
At the present level of development it seems to us that no definitive decision 
can be made as to whether the observed luminosities are being emitted from 
conventional neutron stars or their strange counterparts. 



In summary, we find that really stringent conclusions about the behavior of 
superdense matter on the basis of the observed cooling data cannot be made 
at present. However this situation would change if besides the luminosity and 
the age also the mass of just a single pulsar could be determined! The rapid 
discovery pace of newly observed pulsars in combination with the enormous 
progress that has been made in observational X-ray astronomy are very en- 
couraging in this respect. 
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